function s = sigmasqonxisq

s.sq1 = @sq1;
s.sqonxisq = @sqonxisq;

function y = Ak(k,T)
y = sqrt(2/T)/(1+((k-1)*pi/T)^2);

function y = Bk(k,T)
y = (k-1)*pi/T;

function y = sq1(t,T) % checked
y = sqrt(2)*(1 - (t - 2/exp(T) + 1/(2*exp(2*T)) + 2./exp(t) - 1./(2*exp(2*t)))/T);

function y = sqonxisq(k,t,T) % checked

y = sqrt(2)*(-Ak(k,T)^2*(1/(2*exp(2*T)) - 1./(2*exp(2*t)) + cos(2*Bk(k,T)*T)/2 ...
    - cos(2*Bk(k,T)*t)/2 + (2*cos(Bk(k,T)*t))./exp(t) - T*(Bk(k,T)^2/2 + 1/2) ...
    + t*(Bk(k,T)^2/2 + 1/2) - (2*cos(Bk(k,T)*T))/exp(T) ...
    + (sin(2*Bk(k,T)*T)*(Bk(k,T)^2 - 1))/(4*Bk(k,T)) ...
    - (sin(2*Bk(k,T)*t)*(Bk(k,T)^2 - 1))/(4*Bk(k,T))));
 
 















